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The processes of Z?(*^-niesons and Ac-baryons production in e'''e~-annihilation 
at 10.58 GeV and 91.18 GeV energies are concerned. At the 10.58 GeV energy 

f— ^ I the production of charmed particles via the iJ-mesons decays is also concerned. 
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Scaling violation of the fragmentation functions is calculated at the NLL-accuracy. 
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Nonperturbative fragmentation functions are retrieved from the experimental data 
of S-factories and are approximated by simple analytic expressions. It is proved, 



1^ I that the difference between nonperturbative fragmentation functions of mesons and 

(S i' baryons can be easily explained by quark counting. 
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O ' The production of heavy quarks (charm and bottom) in the high-energy colhsions and 



their subsequent hadronization are important processes in particle physics. Apart from being 
of interest per se, these processes are significant for some practical purposes. For instance, 
a light Standard Model Higgs boson is expected to decay preferably into heavy quark- 
antiquark pairs, which then fragment into heavy hadrons. The direct production of heavy 
quarks would, thereby, be the main background process. Understanding of hadronization 
processes is also important for the precise determination of the mass of the decaying particle. 
Finite heavy-quark mass niQ, providing a natural infrared cutoff, allows to carry out the 
calculations to a large extent using perturbative QCD. Nonetheless, differential momentum 
distributions of the heavy hadrons produced in high-energy collisions are sensitive to the 
large logarithms In s/rriQ, where s is the center-of-mass energy squared. Since for modern 
colliders ^/s ^ mg, these logarithms threaten the convergence of perturbation expansions. 
Fortunately, it is possible to show that up to the power corrections in itiq/s the cross 
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section factorizes into mass-independent hard partonic cross sections convoluted with so- 
called fragmentation functions describing the probability for the parton to fragment into a 
particular heavy hadron [ll, l2|. 

Let us consider the inclusive production of a heavy hadron H via the decay of a vector 
boson \^ = 7*, Z" produced in the e^e~ annihilation: 

e+e- ^V{q)^H{pH)+X. (1) 

It is convenient to introduce the scaling variable x for expressing the heavy hadron energy 
in the center-of-mass frame, 

X = ^ ^^- ^. (2) 

Experimental data are often presented in terms of the scaled momentum Xp = 



\Ph\ /ys/4: — m'jj. But due to large energy of interaction both scaling variables are almost 
indistinguishable. 

As stated above, at leading power in ttiq/s the differential production cross section can 
to all orders of perturbation theory be factorized as 

where daa/dx is the cross section for producing a massless parton a with the scaled energy x 
after subtracting the collinear singularity in the MS factorization scheme, and fragmentation 
function Da/H gives the probability for a parton a to fragment into a heavy hadron H 
carrying a fraction x of the parton's momentum. Da/H also depend on the factorization 
scheme, but the convolution of the two term is not, so that the physical cross section is 
prescription independent [ij. 

The factorization formula ([3]) separates the dependence on the heavy-quark mass rriQ 
from the dependence on the center-of-mass energy which is contained in the partonic cross 
sections. These hard cross sections can be calculated in the massless approximation. All 
dependence on the resulting hadron resides in the fragmentation functions which are process- 
independent non-perturbative quantities. Like the parton distribution functions (PDF's) 
fragmentation functions must be measured at some scale and their values at any other scale 
can be obtained by solving the DGLAP evolution equations 3|-l5|]. 

The fragmentation functions incorporate long-distance, non-perturbative physics of the 
hadronization process in which observed hadrons are formed from the partons. For proper 
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However, such a relation can be expected to hold on 



understanding of hadronic uncertainties one needs to separate short p ^ q and long 
p ~ Aqcd distance effects. The most popular approach is to factorize the fragmentation 
function into perturbative and non-perturbative components: Da/n = D'^Th ® ^cf/H SrlS]. 
The first component is identified with the so-called perturbative fragmentation function, 
D^f^{x^mQ, fi) = Da/Q{xjriQ, n) while for the non-perturbative component a model such 
as Kartvelishvili et al. |9|, [lO| or Peterson et al. ll| is adopted. 

Perturbative fragmentation functions Da/q describe the probabilities for partons a to 
fragment into an on-shell heavy quark Q. They are relevant for the discussion of inclusive 
heavy-quark production, where one sums over all possible hadron states H containing heavy 
quark Q. Quark-hadron duality suggests that 

Y,Da/H{x,mQ,^) = Da/Qix,mQ,iJ,) + o( ^^^ ). (4) 

y if X is not too close to 1, so that 
the scale mQ^l — x) is in the short- distance regime [12|. It is important to mention, that 
such an ansatz spoils the proper factorization of short- and long-distance contributions. 
For X — )■ 1 perturbative fragmentation function itself contains long-distance contributions 
from logarithms of momentum scales of order r?7,Q(l — a;) ~ Aqcd- Such logarithms are 
responsible for soft gluon emission and are not controllable in the perturbation theory. The 
previous attempts to resum the ln"(l — x) terms in the fragmentation functions have thus 
led to unphysical negative values in the a; — )■ 1 region [13']. As a matter of fact, the Landau 
singularity of the perturbative coupling at small transverse momenta leads to a branch- 
point in the resummed expression and produces the negative behavior at 1 — x ~ Aqcd/''^q- 
Soft-gluon resummation therefore suggests that the non-perturbative phenomena become 
dominant when x > Xbr ~ 1 — Aqcd/'^q- The invariant mass of heavy quark and soft gluons 
emitted can be estimated as mg^l + (1 — Xbr)/xbr) ~ ""^g/ll ~ Aqcd/'^q) ~ ttt^q + Aqcd- It 
means that soft gluons revealing themselves at Aqcd scale can play appreciable role in the 
hadronization process. 

To incorporate non-perturbative hadronization effects into the heavy-quark fragmenta- 
tion functions is the main objective of this work. Thus, let us not to resum to all orders of 
perturbation theory the long-distance terms in the perturbative contribution. The motiva- 
tion against resummation was adduced above. The non-perturbative fragmentation function 
will be numerically retrieved from experimental data. Being retrieved from the i?-factories 



data on the D* production at the T(4S') energy, this function will allow to describe the 
ALEPH data at Z-boson peak with the reasonable precision. Apart from testing evolution 
of the D* fragmentation the difference between meson and baryon fragmentation will be 
studied. It will be shown that in the x — )■ 1 region this difference is in agreement with the 
premises of Kartvelishvili et al. model. 

II. THEORETICAL PRELIMINARIES 

A. Perturbative fragmentation function and QCD evolution 



With the use of factorization relation for fragmentation function and the assumption that 

D^f^{x,mQ, fi) = Da/Q{x,mQ, fi) eq. ([3]) can be rewritten as follows: 
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where daq/dx is the heavy quark differential inclusive cross section. 

The MS fragmentation functions Da/Q obey the DGLAP evolution equations 
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where the standard two-loop expression for ^^(yu) is used. Then the perturbative expansion 
for the Altarelli-Parisi splitting functions Pba would have the following form: 



>(o) 



Pba (x, a,(/i)) = a,(/x)Pr(^) + «s(/^)nl^(^) + Oiat), 
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where the P^J are 



<(-) 


= Cp 




1 + x^ 3 

(l-a;)+ + 2 


5{l-x 


p^^H^) 


= 2Ca 


X 1 — X 

(1 — x)+ X 


p>) 


= Up , 
X 


<H-) 


= Tp 


X 


' + {l-xf 


•) 



+ x(l-x) + (^-|^)<5(l-x) 



(9) 



The NLO splitting functions P-j (needed to achieve NLL accuracy) have been computed in 
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18| and are too lengthy to be replicated here. 



The initial conditions for the MS fragmentation functions were first obtained at the NLO 
level in [l|. They are given by 



DQ}Q{x,mQ,fio) = 6{l-x) +as{fio)d''Q\x,mQ,fio) + 0{a1) 



(10) 



(other Da/Q are of order a^), where 
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Although the sum in expression (|5]) runs over all types of partons, -Dg/g is a^-suppressed 
with respect to Dqjq while other Da/q being a^-suppressed. So in the following let us keep 
only direct component Dqjq since it is quite sufficient for the purposes of current work. 
Thus, expanding the convolution, for heavy quark spectrum one obtains 



daq . /- , [^ dz daq 

-{x,y/s,mq,jj) = 



dx 



z dz 



x/z,y/s)DQ/q{z,mq,fi), 



(12) 



Q 



where NLO expression for the partonic cross section from pj| should be used 

-{x,^/s) = (5(1 - x) +as(yu)aQ (x, a/s). 
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The procedure outlined above guarantees that all leading and next-to-leading logarithmic 
terms of quasi-collinear origin (terms of the form {o.s\og{(Pjm^))'^ and as{as\og{q^ /rn^))'^ 
respectively) are correctly resummed in the cross section 1|. 

For the subsequent analysis it is convenient to turn to the Mellin moments of the quantities 
involved. The Mellin transformation f{N) of function f{x) is defined as 

r-l 



f{N) 



dx X f{x) . 



(14) 



In Mellin space the evolution equations ([6]) take the simple form 
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This equation at NLO level was solved analytically in |l| and for the direct component 
Dq/q one has: 

DQ/Q{N,mQ,fi) = E(A^,/i,/io)£'Q/Q(A^,mQ,/io), 



E{N, /i, /io) = exp < In — ^^ ^P-; h 



as(/i) 27r6o 



+ 



asjl^o) - Qs(/^) 
47r2feo 



^■^('V) - ^^('V) 



(16) 



Defining 

JO aa; 

one has the following expression for the NLO distribution: 

a,(iV, v^) = agiN, ^) E{N, fi, /io) D'^)q{N, /xq, mg) 



(17) 
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Both an and dg contain terms proportional to the ^^/(l 
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and 



as [ln(l — 2;)/(l — x)]_^_, associated to the emission of a soft gluon. These terms give 
rise to a large- A^ growth of the corresponding Mellin transforms 
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Leading a" In" A^ and next to leading a^ In" A^ logarithmic contributions were re- 



summed to all orders of perturbation theory in 20|. Opposed to fixed order calculation, 
which leads to finite and positive fragmentation function at almost all values of x except 
the X — > 1 region (where the 5- function from the initial condition becomes apparent), NLL 
resummed result exhibit pathological negative behavior when x approaches 1. The reason 
is that the Landau singularity of the perturbative QCD coupling at small transverse mo- 
menta leads to branch-points in the resummed expression for the initial condition and the 
coefficient function. 



In the initial condition the singularities start at the branch-point 



tL / 1 \ _ , /^O 
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while for the coefficient functions the branch-point is 
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< = exp ( 7-^ ) ^ T^ • (21) 

The previous attempts to restore the physical behavior of fragmentation functions con- 
sisted in introducing a tower of power corrections to A^ represented by the replacement 

N ^N ^ ^ ^^^Z , (22) 



'ini 



in the initial condition and 

in the coefficient function. It is easy to see that / being more or equal to 1 unphysical region 
is unreachable. But it is important to keep in mind that there is no rigorous justification 
for the replacements ( l22l) . ( l23l) . 

B. Non-perturbative fragmentation function 

The most popular parameterizations for the non-perturbative fragmentation functions 
are Peterson et al. ll| and Kartvelishvili et al. (KLP) [9|, [lO|. The former has a form 
of the heavy quark propagator and does not depend on the hadron produced. Thus let 
us concentrate on the latter one. It is based on the Gribov-Lipatov "reciprocity relation" 
between D^,^ and the distribution function of quark Q in hadron H [^ : 

^Q%(^) '= /^(^), (24) 

where z = Ph/Pc is the hadron momentum fraction with the respect to the heavy quark mo- 
mentum. The expression for /^, was found out in [26] on the basis of Kuti-Weisskopf model. 
The parametrization obtained is significantly related to the Regge trajectory parameters of 
the QQ-system and has the following form: 

/£.(-) = rn""'?;n?"°'\ ^-°(i " -)'"-• (25) 

r(l -aQ)r(l + 7a/ - ag) 
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where aq = 1/2 is the intercept of the hght quarks trajectory p, u, f, A2, oq is the intercept 
of the leading trajectory for the QQ-system and jh is a parameter determining the behavior 
of the distribution function for a; — > 1. Its value originates from the q^~'^^ diminution of the 
form- factor as it is known that q^^'^^'> diminution of form-factor leads to the (1 — x)^A; — 1 
behavior of the distribution function at large x. Assuming the universality of the sea quarks 
distribution in all mesons one gets 7*/ = 3/2. In much the same way for A^-baryons the 
following expression was found out: 

■'^'^ ' r(l-aQ)r(l + 7B-2a,) "- ' ' ^ ^ 

where 7b = 3 and the factor 2 is related to the number of valent light quarks in baryon. 

Gribov-Lipatov "reciprocity" relation with known distribution functions (125|) . (126|) deter- 
mine the behaviour of the fragmentation functions in the x — )■ 1 limit. For the small values 
of X the condition for the fragmentational approach to be valid (p ^ ttiq) is not satisfied. 
Nonetheless, if at small x values the production of heavy hadrons mainly depend on their 
wave functions, then coinciding x — > 1 asymptotes of fl25l) and fl26l) prognosticate similar 
behaviour of the momentum distributions of mesons and baryons in this region. 

There is still some uncertainty in the value of a^- Theoretical investigations [27| based 
on Regge trajectory systematics result in the value for ac in the range between —2.0 and 



—3.5. It is slightly more than the value of ac ~ —3 -. — 4, obtained in 



with the use of 



the QCD sum rules. Another way to determine Uc by the value of the heavy quarkonia wave 



function in the center point leads to the value —3.5 ± 0.6 |29 |. 



III. CHARM HADRONS DATA FITS NEAR THE T(45) 

High quality data on the charmed hadron production is provided by BELLE, BABAR and 



CLEO collaborations 



21I-I23 



23|. One thus has the opportunity to perform a more accurate 
fit for the non-perturbative initial conditions. Furthermore, it gives us the possibility to 
test the evolution of the fragmentation function from the center-of-mass energies of 10.6 to 
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91.2 GeV, using charm data from the LEP experiments 

For the perturbative component of fragmentation function we use expression f fTSj) with the 
NLO initial condition ( ITOj) . the NLO partonic cross section (fT3|l and the NLL evolution ( !T6|) 



33|. As stated above we do not perform NLL Sudakov resummation for the initial conditions 
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and coefficient functions retaining these long-distance contributions for the non-perturbative 
component. 

Several parameters enter the perturbative calculations. First of all these are the initial no 
and the final fi evolution scales which allow variation by a factor of order of 2 around rric and 
^/s respectively. We set them to be /io = 2mc and /i = ■\/s/2 as this values allow to perform 
the most successful fiting. The center-of-mass energy ^i is equal to mx(45) = 10.58 GeV. 
We shall use the pole mass for charm quark and fix it to be rric = 1.6 GeV. As the 6-quark 
mass lies between fiQ and /i it creates a threshold on the different sides of which different 
number of active fiavours enter the evolution operator. We set nib = 5.0 GeV. Experimental 
value of as{mz) = 0.119 points to the value of Aqcd equal 0.226 GeV. 

There are several ways to retrieve the non-perturbative component. DGLAP equations 
allow to perform the evolution from a larger scale to a lower one as well as in the opposite 
direction. But then one needs to extract the non-perturbative function from its convolution 
with the partonic cross section and the perturbative initial condition. This procedure is 
connected with the inverse Mellin transform of the moments of experimental data divided 
by the moments of the the partonic cross section and the perturbative initial condition. Such 
calculation performed numerically faces problems with the integral's convergence developing 
into unphysical negative values of the non-perturbative function in the low x region. It was 
tested that resulting non-perturbative function obtained by such way does not permit to 
successively reproduce the experimental data by the reverse procedure. 

Anther method, proposed in this work, is to represent the function required as a linear 
combination of functions which have a simple analytic form of Mellin transform. Further, 
expansion coefficients can be determined by fitting to experimental data in the Mellin space 
as well as in the a;-space. The most general choice is to retrieve the non-perturbative 
fragmentation function bin-by-bin, choosing the number of bins n to be not more than the 
number of experimental points to avoid overdetermined system for the coefficients. So, let 
us define 

5-W^|:e,e(.-ii^)e(l-.), (27) 

The corresponding Mellin transform is 

The approximate number of points in the BELLE and BABAR data is 50, in the CLEO 
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FIG. 1: Nonperturbative fragmentation func- 
tions of the _D* -mesons extracted from several 
datasets. 



FIG. 2: Averaged nonperturbative fragmenta- 
tion function of the Z3*-mesons. 



data — 20. Thus to take 20 bins for the non-perturbative function required seems to be a 
reasonable choice. 

Fit to the experimental data in the A^-space can be performed as well as in the x-space. 
The Cj coefficients obtained by both techniques are in a good agreement. There are 4 data- 
sets for D* production at T(4S') energies. Two of them presented by BELLE and two by 
CLEO, they regard to the D*^ and D*^ production. The non-perturbative fragmentation 
functions extracted from them are plotted in Fig. [1] As we assume that there is no difference 
between the D*^ and D*^ fragmentation these functions can be averaged to get the D* non- 
perturbative fragmentation function plotted in Fig. |2l The weights in this average were 
selected proportional to the statistics gathered for each data-set. 

The same procedure was carried out for the A^-baryon production at BABAR and 
BELLE. The corresponding plots for the non-perturbative fragmentation functions are pre- 
sented in Fig. [3]and|ll 

The mesonic and baryonic fragmentation functions obtained in such way do not reveal 
significantly different behaviour aX z "^ 0.5 (Fig. [5]). Furthermore, best fits by c-x~° function 
in the x < 0.5 region for both cases result in close values of a. These values being equal to 
—3.7 and —3.8 are in a pretty good agreement with the previously mentioned predictions 
for Oc- What concerns the 2; — ?■ 1 behaviour, the difference in it is in agreement with the 
Gribov-Lipatov "reciprocity" relation. 

For further phenomenological analysis it is convenient to find some simple expression 
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FIG. 3: Nonperturbative fragmentation func- 
tions of the Ac-baryons extracted from several 



datasets. 



approximating the numeric data obtained: 
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FIG. 4: Averaged nonperturbative fragmenta- 
tion function of the Ac-baryons. 
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20.1z''-'(l - z) + 2.7710V^(1 - z)\ 
D^^iz) = 72.9z^-'^(l-zf + 2.93 lO^z^^ (1- z)^ + 10^ z^'^n-zf. 



(29) 
(30) 



Each of them has the same z — )■ and 2; — )■ 1 asymptotes as the corresponding KLP function. 



IV. CHARMED HADRON PRODUCTION IN ANNIHILATION PROCESSES 

According to the formula ([5]) moments of the momentum spectrum of the particles pro- 
duced are equal to the product of moments of the differential partonic cross-section flTSj) and 
moments of the corresponding non-perturbative fragmentation function. To return to the x 
variable the inverse Mellin-transform should be performed by integrating over the vertical 
line in a complex plane: 
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(31) 



As in a fixed-order calculation the Landau pole does not appear it is possible to use any pos- 
itive value of 7. Coincident results obtained at different values of 7 prove the independence 
on its value. 

Non-perturbative functions (l29l) and (130|) allow to reproduce the experimental data from 
S-factories with good precision, see Fig. [6l [7] and [H 
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FIG. 5: Nonperturbative fragmentation func- 
tions of _D*-mesons and A^-baryons together 
with the approximative expressions (I29|) and 




FIG. 6: _D*+-meson momentum distribution 
from the c-quark fragmentation at ^/s = 
10.58 GeV energy together with the Belle and 
CLEO experimental data. 

(T,Q, dcr/dx 




FIG. 7: £'**'-meson momentum distribution 
from the c-quark fragmentation at ^/s = 
10.58 GeV energy together with the Belle and 
CLEO experimental data. 



FIG. 8: Ac-baryon momentum distribution 
from the c-quark fragmentation at ^/s = 
10.58 GeV energy together with the BaBar and 
Belle experimental data. 



Evolution to the scale \/s/2 = 45.6 GeV and non-perturbative expression fl29|) are used 
to obtain the momentum distribution of i5*-mesons at the Z-boson peak. The predicted 
spectrum together with the ALEPH and OPAL data is presented in Fig. [91 These distribu- 
tions coincide with the reasonable precision. As the only difference in calculations for the 
10 and 90 GeV energies was in the final evolution scale, one can state that factorization 
relation is valid in this energy range. 
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FIG. 9: £'*+-nieson momentuni distribution 
from the c-quark fragmentation at ^/s = 
91.18 GeV energy together with the ALEPH and 
OPAL experimental data. 



FIG. 10: Predicted Ac-baryon momentum dis- 
tribution from the c-quark fragmentation at 
yi = 91.18 GeV energy. 



For the Ac production at 91.2 GeV energy the same evolution in the perturbative com- 
ponent is used. The non-perturbative effects are described by the expression ( 130|) . Unfortu- 
nately there is no experimental data on the Ac production at Z-boson peak. Our prediction 
for the Ac momentum distribution is presented in Fig. [101 

A considerable part of Z^'s is produced indirectly through D* decays. The simple spin- 
state counting estimation (2 J -|- 1) leads to the factor 3 enchantment for D* production. 
The experimental value obtained at Z-boson peak amounts to only 1.4 3l|- We will use 
experimental value and assume that D*^, D*~ , D*^ and Z)*° 3^ are produced with equal 
probabilities. 



Following the approach of 30|] we assume that the D meson non-perturbative fragmenta- 
tion function is the sum of a direct component, which is isospin invariant plus the component 
arising from the D* decay. 

The decay D* — t- Dtt takes place very close to the threshold. The momentum of D-meson 
in the D*-meson rest system 
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(32) 



is sufficiently small to be neglected. Thus D has the same velocity as the D*, and their 
momenta are thus proportional to their masses. So the component of the D-meson fragmen- 
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tation function arising from D* — ^ Dn decay is given by 



z^-(.) = z^f f.^Vfl-.^)^, 



where D^ {z) is the non-perturbative fragmentation function of D*-meson fl29|) . The integral 
of the expression (l33ll equals 1, so if should enter the D-meson fragmentation function with 
the weight proportional to the D* — t- Dvr decay probability and probability of the D* 
production. 

What concerns the D* — )■ D^ decay, momentum of the D in the D* frame is non- 
neghgible: 



P 



2m£)* 



135 . 



(34) 



The D momentum in the laboratory frame is given by a Lorentz boost 



P = 'J {p' cos 9 + f3e'), 



(35) 



where (3 is the velocity of the D*-meson, 7 = l/-\/l — [3'^, e' = {mj^, + m|,)/(2m£i.) — 
energy of the D-meson in the D* rest frame and 6 — its decay angle with respect to the 
D* direction. Denoting momentum and energy of the D* in laboratory frame by p* and e* 
respectively one obtains 
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Introducing variables 



P 



P 



^D 



p* _ p* 



Pmax \/'S/4 — m?j-), 

contribution of the D* — ^ D7 decay to the D production can be written down as 



J-l 2 \ Pmax 



As in the previous case the integral o 
The branching ratios involved are 



this expression is normalized to unity. 



3^: 



(36) 



(37) 



(38) 



BrD*+^DOn+ = 67.7 ±0.5, 
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Sr^.+^m^o = 30.7 ±0.5,% 

BrD*+^D+-y = 1.6 ±0.4,% (39) 

Btd^o^do^o = 61.9 ±2.9,% 

Bro^o^DOy = 38.1 ±2.9.% 

Finally non-perturbative fragmentation functions of D-mesons can be written down as 
follows: 

Dj^^iz) = n^^(Df(^)±c[5r^*+_^+,D^^(z)± 

± BrD*+^D+noD''^iz)]) (40) 

and 

D^\z) = n'^'iD^iz) + c[Bro*o^no^D''^iz) + 

+ {BrD*+^DOn++BrD*o^Do^o)D^''{z)]), (41) 

where c = 1.4 — ratio of probabilities to fragment into D* and D-mesons, coefficients n^ 
and n^ provide normalization of fragmentation functions to unity: 



n 



D+ 



DO 



^1 ± c^Brn'O^oo-r ± 5r£,.+^D0^+ ± Bro^o^o'^no))'^- (42) 



n 

In order to obtain momentum spectra of D-mesons expression (I3T1) is used. The result 
for i5+-mesons is presented in Fig. [TTl for D^-mesons — in Fig. [121 Both distributions are 
in a good agreement with experimental data. The D^ spectrum is slightly softer then the 
D^ one because of the larger probability of the D* — t- D^X decay. 

V. CHARMED HADRON PRODUCTION IN S-MESON DECAYS 

Let us now consider charmed particles production in the 6-quark decays. The energy of 
10.58 GeV corresponds to the T(4S') resonance which decays into a i?i?-pair with almost 
unitary probability. S-mesons from such decays are nearly at rest as their mass niB = 
5.28 GeV ^ -\/s/2. Neglecting 6-quark motion within the meson the c-quark spectrum from 
B decay can be easily found. The charmed hadron spectrum can then be written down as 

'""w^ r^ fowler f-). (43) 



dx Jx z \ dz I '^ \z 
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(Tfgfdcr/dx 



cr^Jfdcr/dx 





FIG. 11: Z3+-meson momentuni distribu- 
tion from the c-quark fragmentation at y/s = 
10.58 GeV energy together with the Behe and 
CLEO experimental data. 



FIG. 12: Z?'^-meson momentum distribution 
from the c-quark fragmentation at -y/s — 
10.58 GeV energy together with the BeUe and 
CLEO experimental data. 



where dab-^c/dz is the c-quark spectrum obtained by analysis of weak decays h — )■ c + Ivi, 
b ^}- c + qq, b ^ c + cs and b ^ c + cs. 

We neglect here any perturbative fragmentation function as at such low energy it should 
not sufficiently differ from the ^-function. The non-perturbative component of D* production 
is described by the expression (!29ll . 

In experimental data value x = 1 corresponds to the largest possible momentum of the 
D*-meson prodused at 10.58 GeV energy: 



P: 



D* 



s/A - m%, = 4.88 GeV 



(44) 



while the largest possible momentum of the D*-meson from the B decay is equal to the 
largest c-quark momentum 



/^niax 



\ 



Ami 



-ml = 2.24 GeV. 



(45) 



This value actually coincides with those calculated via the hadron (not parton) masses: 



rvaax. 



\ 



2 \2 



(m| + m \) 
4m| 



VTi 



D 



2.26 GeV. 



(46) 



Thus the momentum distribution of Z^-mesons from the B decays is located in the region 
^ < Pmax/Pmax = 0.46. The distribution concerned and the experimental data are plotted in 
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o",„' dcr/dx 





FIG. 13: _D*-meson momentuin distribution 
from the 6-quark fragmentation at ^/s = 
10.58 GeV energy together with the Belle ex- 
perimental data. 



FIG. 14: Z3*-nieson momentum distribution 
from the 6-quark fragmentation and recombi- 
nation with the light quark from _B-meson at 
V^ = 10.58 GeV energy together with the Belle 
experimental data. 



Fig. [131 In this only plot the normalization of the experimental data was changed for points 
in the low x region to be in agreement with the predictions. In the high x region there is 
some extra contribution not predicted by the fragmentation model. We suppose that this 
contribution arises from the interaction of the charm quark with the light valent quark from 
the 5-meson. As this quark has larger momentum then a sea quark the process concerned 
contributes to the high x region. We assume that the momentum distribution of D*-mesons 
produced in this process coincides with the c-quark distribution. This assumption allows to 
obtain better agreement with the experimental data (Fig. [T^ . 

Let us now proceed to the Ac-baryon production. To obtain Ac-baryon momentum 
distribution one should convolute c-quark spectrum with the non-perturbative fragmentation 
function fl30l) . The largest possible momentum of the A^-baryon prodused at 10.58 GeV 
energy is equal 



Pmax = \/s/4 - ml^ 



4.76 . 



(47) 



When finding largest A^-baryon momentum from the B decay it is important to mention 
that baryon production is always accompanied by the any-baryon emergence. The lightest 
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CTigfdcr/dx 




FIG. 15: Ap-baryon niomentum distribution from the 6-quark fragmentation at ^/s = 10.58 GeV energy 
together with the BaBar and Belle experimental data. 



of them is anti-proton. Thus the largest possible A^-baryon momentum equals to 



/^max 



\ 



4m| 



"^A. 



2.02 . 



(48) 



max' 



This value differs significantly from the p^ 

Thereby the distribution sought for occupies the region x < Pmai^VPmax — 0.42. It 
reasonably coincides with the experimental data (Fig. [T5|) . As opposite to the meson 
production sea wrf-diquark capturing in needed for the Ac formation is needed. Tat is why 
valent quark from the 5-meson does not change the prediction of fragmentation approach 
in case of baryon production. 



VI. CONCLUSION 



In this article charmed hadron production in the vast range of energies was concerned. 
During this the energy dependence was contained only in the perturbative component of 
fragmentation function, for which the NLO-expression was used. Good agreement with the 
experimental data points to the fulfilment of the factorization assumption in the 10 to 90 
GeV energy range. 

It is important to mention that non-perturbative fragmentation functions were one and 
the same for all energies except D*-mesons production in B decays. It means that separation 
of non-perturbative phenomena was carried out correctly. Indeed the perturbative part is 
relevant for charmed quark production only, the non-perturbative — only for its transaction 
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to the final particle. 

The difference in meson and baryon production in the x — )■ 1 region dues to the different 
non-perturbative fragmentation functions. The difference in them by-turn is explained by 
quark-counting. This phenomenon takes place in large x region only, where non-perturbative 
phenomena is most significant. 

In the X — ;■ region non-perturbative fragmentation functions of charmed mesons and 
baryons coincide. This is relevant to the one and the same behavior of charm quark dis- 
tribution functions in mesons and baryons in the low-x limit. Moreover, the parameter a 
which determines the low-x behavior is close to the estimations of the ac — interception of 
the charmed Regge trajectory. 

At the low energy scale which corresponds to the S-meson mass a discrepancy in the D* 
momentum distribution was found. The origin of it is thought to due mainly to the influence 
of a valent quark from decaying B-meson. Such discrepancy does not occur in Ac production 
which is in a pretty good agreement with the experimental data. The rri/i/i-corrections to 
the factorization relation also worth mentioning at the rriB energy scale. 

Author would like to thank Prof. Likhoded A.K. and Luchinsky A.V. for useful ideas and 
discussions. The work was financially supported by Russian Foundation for Basic Research 
(grant#10-02-00061-a) and grants of the president of Russian Federation MK-140.2009.2 
and MK-406.2010.2. 
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